Driving-dependent damping of Rabi oscillations in two-level semiconductor systems 
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We propose a mechanism to explain the nature of the damping of Rabi osciUations with increas- 
ing driving-pulse area in localized semiconductor systems, and have suggested a general approach 
which describes a coherently driven two-level system interacting with a dephasing reservoir. Present 
calculations show that the non-Markovian character of the reservoir leads to the dependence of the 
dephasing rate on the driving-field intensity, as observed experimentally. Moreover, we have shown 
that the damping of Rabi oscillations might occur as a result of different dephasing mechanisms for 
both stationary and non-stationary effects due to coupling to the environment. Present calculated 
results are found in quite good agreement with available experimental measurements. 
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Localized semiconductor systems exhibiting few dis- 
crete energy levels ("artificial atoms"), such as specially 
selected donor impurities and quantum dots (QDs), are 
prospective candidates to play the role of basic build- 
ing blocks for quantum information processing. In par- 
ticular, a two-level semiconductor system may exhibit 
Rabi oscillations (ROs) of its population when coupled 
to a drivingfield, so that it may be coherently controlled 
[H, [3, H, There are a number of dephasing mecha- 

nisms for localized semiconductor systems, some of which 
arc essentially non-Markovian so that one needs to take 
into account memory effects as well as the back-action 
of a dissipative reservoir on the radiating system. For 
example, a dephasing caused by spin-spin coupling be- 
tween neighboring QDs or carriers captured in traps in 
the vicinity of a QD was shown to lead to non-Markovian 
dynamics @, Such reservoirs have correlation times 
comparable with the typical decoherence time of the de- 
phasing system. Also, the dephasing due to coupling 
with phonons was shown to lead to non-Markovian fea- 
tures in the dynamics of a two- level systems (TLS) Q. 
Carriers and excitons in localized semiconductor systems 
may be coupled not only to localized neighboring states, 
but also to delocalized ones This diversity of dis- 
sipation channels has led to a number of novel features 
in such systems' dynamics. In the present work we fo- 
cus our attention on one peculiar phenomenon which has 
caused and is still causing much controversy, namely, the 
damping of ROs due to the increase of the driving-pulse 
area which is an observed feature of coherently excited 
localized semiconductor systems [l|, 0, [4, 5] . A number 
of mutually contradicting explanations was suggested for 
it. One of these is that such a dephasing is due to the 
system's interaction with a non-Markovian reservoir of 
phonons 0. However, the dephasing process takes place 
even when the coupling with phonons is negligible 
Driving-dependent damping of ROs was proposed to oc- 
cur as a consequence of excitations of bi-excitons in the 
QD [10], although damped ROs are also observed when 



there is no possibility for the bi-exciton excitation 
Recently, it was demonstrated that the experimentally 
observed 0| intensity-dependent damping of ROs can be 
reproduced by introducing into the standard Bloch equa- 
tions a dephasing rate dependent on the driving-field in- 
tensity [llj . On the other hand, although there is an ex- 
perimental confirmation of a driving dependence of the 
dephasing rate an intensity- independent dephasing 
rate has also been measured Q ■ 

Based on this controverted scenario, in the present 
work we propose to shed some light on this matter 
by studying a simple TLS excited by a classical coher- 
ent field and coupled to a general dephasing reservoir. 
Within a quite general and straightforward approach, we 
demonstrate that a driving-field dependent damping of 
ROs stems from various relaxation mechanisms entering 
into play in different experimental situations. Further- 
more, we show that driving-dependent damping may oc- 
cur whether the reservoir is influenced or not by the driv- 
ing field. To keep it simple, and to focus only on features 
which give rise to the phenomenon in question, we assume 
no population damping of the TLS. This also corresponds 
to the real experimental situation with driving by short 
laser pulses, so that the population damping is negligibly 
slow on the time-scale of the system's dynamics In 
the frame rotating with the driving-field frequency, lol, 
working within the interaction picture with respect to 
the reservoir variables and using the rotating-wave ap- 
proximation (RWA), we describe our problem with the 
following standard effective Hamiltonian, 

Htot{t) = Ho{t) + ha+a-YL{t), (1) 

where the undamped system's Hamiltonian is given by 

Ho{t) = hAa+a- + h[n{t)c7+ + n*{t)a-]. (2) 

Here cr* = |i)(T| are the system's raising and lowering 
operators, the kets |±) correspond to the excited and 
ground states of the TLS, respectively, A = loq — lol 
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is the detuning of the driving-laser frequency u>l from 
the resonance frequency of the TLS transition, with 
the possible addition of a frequency-shift term due to the 
interaction with the dephasing reservoir. The reservoir is 
described by the operator R(t) , which might also depend 
on classical stochastic variables (describing, for example, 
different realizations of the reservoir in each run of an 
experiment ^]), whereas fl{t) describes the shape of the 
driving pulse. 

Let us now assume that the reservoir correlation func- 
tion (R(i)R(T)) satisfies the following general require- 
ments: (R(t)R(r)) ^ 0, when t,T ^ oo, and \t — 
t| ^ oo. If the coupling of the reservoir to the TLS 
is weak, and the reservoir correlation function decays 
with t,T —^ oo, and also with |t — t| — > oo much 
faster than the typical time-scale of the system's evo- 
lution, it is possible to obtain a time-local master equa- 
tion V, for the problem described by the Hamilto- 
nian I!])-©. Following the approach developed in Refs. 
0, we introduce dressed operators describing the inter- 
action with a classical field, i.e., S^(t) = UHt)a^U(t), 
where the unitary " dressing" transformation [15| is given 

by: U{t) ^ ¥exp |-|- /i/o(T)dT|, and ^ denotes 

the time-ordering operator. One may use the time- 
convolutionless projection operator technique or cumu- 
lant's expansion and the Born approximation for the 
"dressed" density-matrix master equation, and then go- 
ing back to the "bare" basis, one obtains the following 
set of Bloch equations with time-dependent coefficients 
iH: 



dt 

dt 



= i[n{t)p^+-n*{t)p+^], (3) 
= {i[^ + {^m-<^}p+- (4) 

where p++ = p±^ = (±|p|=F), p is the density 

matrix of the TLS in the "bare" basis, and the time- 
dependent dephasing rate K[t) and the generalized Rabi 
frequency Q,{t) are defined as 



K{t) = J dT(R(T)R(t)>Z?++(r-i), 



(5) 



to 



n{t) = n{t) ~ J dT(R{T)Ii{t))D+_{T -t) , (6) 

to 

where (t) and D^+{t) are dressing functions In 

the case of a rectangular pulse = for the pulse 

duration, which we will use in further discussions here), 
one obtains 



D++{t) = 



in 



(7) 



where c = A/flji, s = fi/fJ^j, and 17^ = + fi^ is the 
effective Rabi frequency. 

Let us now consider the simplest situation in which 
the driving field interacts only with the localized system. 
In the Markovian limit one has (R(r)R(i)) ~ S{t — t) 
and, therefore, as follows from Eqs. (O-®, one recovers 
the standard system of Bloch equations for a driven TLS 
in the presence of dephasing effects. In this case, the 
dephasing rate, K{t) = k, is constant and independent 
of the driving-field intensity. Then, as expected, ROs 
persist for all values of the field's intensity [cf. dotted 
curve in Fig. 1(a)]. 

For a general non-Markovian reservoir one may write 
the reservoir's correlation function as a sum of a sta- 
tionary contribution K{t — t) and a non-stationary one 
P{T,t) which tends to zero for t = t as t,T ^ oo, i.e., 
(R(T)R(t)) = K{T-t)+P{T,t), where P(r,t) is respon- 
sible for non-Markovian effects at the initial stage of the 
system's dynamics. For the moment, let us ignore effects 
of P{T,t), and consider the Fourier-transform K{w) of 
K{t) = f dwK{w)e-'^'^-'^^^^. FromEq. P, one obtains 



dT / dwKiw)e-'^'"-'^'-^^^-'^D++{T-t) (9) 



to 



for the dephasing rate. Notice that the Markovian ap- 
proximation holds whenever K{ll!) is smooth in the vicin- 
ity of both the frequency ujl and TLS transition fre- 
quency. Moreover, Eq. © indicates that a sufficient 
intense driving-field probes K(lu) away from the ujl fre- 
quency. The K{uj) spectrum may be smooth enough in 
the vicinity of all components of the triplet lul, ujl ± i^/j 
to justify a Markovian approximation for each of them 
Taking into consideration that K{u!) has different 
values at these frequencies, even a Markovian approxima- 
tion for each component of the triplet should yield to an 
intensity-dependent dephasing rate. Therefore, by per- 
forming the Markovian approximation for the different 
components of the triplet in a standard way, for a rect- 
angular driving pulse, one finds from Eq. ([9]) the time- 
independent dephasing rate @: k « ^(c^ -I- \)K{ijJl) + 
js'^[K{ujL + nR)+K{ujL-nR)\. Moreover, when differ- 
ences in values of K{u!) at frequencies ujl, i^l ± ^^_r are 
much smaller than the value oi K{ijJl), one may expand 
K[uj) in the vicinity of wl and obtain 



k^^K{ul) + ——K{uj) 



(10) 



D+-{t) = —{c[l-cos{nRt)]+ism{nRt)}, (8) 



as an intensity-dependent dephasing rate. Here we notice 
that Brandi et al [Tlj have used an intensity-dependent 
recombination rate as in Eq. (jlOp to model experimental 
measurements on ROs in a QD semiconductor TLS, and 
found good agreement with the excitonic photocurrent 
data as measured by Zrenner et al Also, from Eq. 
(pl), one may use the same approximation as before in 
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obtaining Eq. (jlOp . and find 

for the generalized time-independent Rabi frequency. 

Now we apply the developed approach in order to 
obtain a quantitative understanding of the experimen- 
tal measurements by Zrenner et al 2] and Wang et al 
3. Figure 1 displays the present results correspond- 
ing to the solution of the Bloch equations with the 
driving-dependent dephasing rate and generalized Rabi 
frequency [see Eqs. pil)) and (fTTj) ] chosen in order to give 
the appropriate ROs as found in the experimental mea- 
surements 0, 01 • One clearly notes the excellent agree- 
ment with the excitonic photocurrent measurements of 
Zrenner et al [Fig. 1(a)] and photoluminescence mea- 
surements by Wang et al ^ [Fig. 1(b)]. One needs to em- 
phasize that, with respect to the effects stemming from 
the stationary part of the reservoir's correlation function, 
the particular form of the K{uj) function is of no impor- 
tance as long as it satisfies quite general requirements 
as mentioned before. In the present calculation only the 
value of the K(uj) function at the point ujl and two of 
its derivatives are of importance [cf. Eqs. pH]) and pT|) ]. 
These are the only "free" parameters to match the ex- 
periment. Moreover, apart from the value K{u>l), only 
the second derivative of K{u)) at the point uj^ plays a sig- 
nificant role, and we have actually used essentially this 
parameter to match the experimental data. 

We now consider that the coherent driving-pulse ap- 
plied to the TLS may also influence its surroundings. If 
the action of the driving field on the system surround- 
ings is weak, the K(t — t) stationary contribution to the 
reservoir will essentially have the same dependence on the 
driving-field intensity as described above. We note that 
the driving-pulse action on the reservoir may also give 
rise to non-Markovian effects stemming from the P{t, t) 
non-stationary part of the reservoir's correlation func- 
tion, and that observable manifestations of these effects 
may be very similar to those described above. Let us 
illustrate it with a simple model of a bosonic reservoir 
driven by the same rectangular pulse that is applied on 
the TLS under investigation. Using the RWA, one may 
describe the whole "TLS -I- reservoir" system with the 
following Hamiltonian 

Hiit) = Ho{t) + Hresit) 

+ ha+a-J2[gjb+e'^^'-'-*°^+h.c.], (12) 
3 

where Hresit) is the reservoir Hamiltonian, 

3 3 

and the gj are interaction constants, the are detun- 
ings of the reservoir modes from the driving field, and 
the rij are Rabi frequencies for every particular reservoir 




Excitation amplitude (arb. units) 



FIG. 1: (a) Rabi oscillations of the photocurrent, at reso- 
nance, as a function of the excitation amplitude. The dot- 
ted line is the solution given by the Markovian Bloch equa- 
tions with the dephasing rate independent of the driving field, 
whereas the solid line corresponds to the solution of the Bloch 
equations with the driving-dependent dephasing rate and gen- 
eralized Rabi frequency given by Eqs. (|10|) and (lllf) . Full 
squares represent experimental data from Zrenner et al , for 
a pulse width of about 1 ps. Here, a vr-pulse corresponds to the 
unit of the excitation amplitude; (b) ROs in the photolumines- 
cence (PL) intensity, at resonance, with full theoretical curves 
corresponding in descending order to pulse widths of 9.3 ps, 
7.0 ps, and 5.4 ps, respectively. Calculations are performed 
with the driving-dependent dephasing rate and generalized 
Rabi frequency as in (a). Full symbols are the corresponding 
experimental data from Wang et al Qj. 

mode (we assume them to be constant, ^j{t) = fij, for 
the pulse duration). Using the interaction picture with 
respect to Hres (t) , one recovers from Eq. the Hamil- 
tonian of Eq. ID) with the following reservoir operator 

= E 93 h + ^] e'(f2^+A,)(to-t) + ^ (14) 

3 

and with the system's detuning shifted due to the in- 
teraction with the excited reservoir, i.e., A — > A — 
2^gj^j / l^j. For the reservoir initially at the vacuum 

3 

state, one obtains 




and the reservoir correlation function (R(T)R(t)) as the 
sum of a stationary part K{t — t) = ^ ^2gi(at+Aj)(t-r) 

3 

with a non-stationary part P{T,t) = (R(T))(R(t)). The 
stationary part K{t — t) produces effects already de- 
scribed above, and here we assume that K(w) is wide and 
smooth enough so that the stationary dephasing rate, Ks , 
is independent of the intensity of the driving field. Then, 
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from Eq. ([5]), one has 



Hit)) dT{TL{T))D++{T-t). (16) 



to 



Note that the non-stationary part of the dephasing rate 
K{t) decays with time [see Fig.[5l^a)], since (R(i)) for 
t — > cx), and that the function (R(i)) may decay slower 
than the stationary part, K{t), of the reservoir's correla- 
tion function as the driving field excites different modes 
of the reservoir in a different way, and the spectral den- 
sity of the reservoir's excitation may therefore be much 
narrower than K{w). Moreover, in experiments on ROs 
in localized semiconductor systems one deals with short 
driving pulses, so that the non-stationary part of the de- 
phasing rate may play a significant part in the system's 
dynamics. Even if one assumes ~ 1 for the time- 

interval of interest, the non-stationary part of the dephas- 
ing rate will be dependent on the driving-field intensity. 
This is a purely non-Markovian dynamical effect produc- 
ing an intensity-dependent damping of ROs [cf. Fig.[2lb)] 
quite similar to those described before. 

The decrease of the dephasing rate with time may be 
responsible for the constant value of the dephasing rate 
as measured after the application of the driving pulse 
in the experimental measurements by Patton et al Q 
[this situation is illustrated in Fig. [Ha)]- Also, it may 
explain the decreased dephasing rate after the applica- 
tion of the pulse as in the experiment by Wang et al Q . 
To conclude, a driving-dependent damping of ROs due to 
the non-stationary contribution of the reservoir's correla- 
tion function may take place for quite general reservoirs. 
Indeed, the nature of the reservoir influences only the 
particular form of P(r, and not its general properties, 
which determine the effect in question. 

Summing up, we have demonstrated that the damp- 
ing of ROs with the driving-field intensity in localized 
semiconductor systems (QDs, shallow donors, etc) is an 
effect of a very general nature, and a consequence of non- 
Markovian effects due to the coupling of the system to a 
reservoir. The exact nature of a reservoir (being an en- 
semble of phonons, other localized systems, traps, free 
carriers in a wetting layer, coupling to bi-excitons or 



higher decaying levels, etc, or a combination of mech- 
anisms) is not of particular importance for the manifes- 
tation of the effect. Similar damping of ROs may occur 
as a consequence of different physical mechanisms. The 
first one stems from stationary properties of the reservoir 
whereas the second one is a purely non-stationary effect 
occurring when the driving field excites the reservoir with 
a decay time of the non-stationary part of the reservoir's 
correlation function comparable to the driving-field pulse 
length. 

The authors gratefully acknowledge partial financial 
support by EU under EQUIND project of 6FP IST- 
034368 and INTAS, and by Brazihan Agencies CNPq, 
FAPESP, Rede Nacional de Materials Nanoestrutura- 
dos/CNPq, MOT - Millenium Institute for Quantum In- 



T3 

to 



(a) 



3 4 
Time (ps) 




2 4 6 

Excitation amplitude (arb. units) 



FIG. 2: Examples of (a) dephasing rate nit) [cf. Eq. 
(|16p ]. and (b) upper-state population dynamics, for the fixed 
time moment corresponding to the end of the rectangular 
pulse, versus excitation amplitude, for a model (R-(i)) ~ 
7^6"^'*"*"-', with 7 = 2ps~^. Dotted lines correspond 
to D++(r — t) = 1 whereas solid lines correspond to the 
D++{t - t) defined in Eq. 0. 
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